The Ehrenfest system and the rest point spectrum for a 
Hartree-type Equation 



V. V. Belov^ M. F. Kondratieva*"^ A. Yu. Trifonov*^ 

Department of Applied Mathematics, 
Moscow State Institute of Electronics and Mathematics, 
Trehsvjatitelsky per., 3/12, Moscow 109028, Russia 

* Department of Mathematics and Statistics, 
Memorial University of Newfoundland, 
St. John's AlC 5S7, Canada 

* Laboratory of Mathematical Physics, 
Mathematical Physics Department, 
Tomsk Polytechnical University, 
Lenin ave. 30, Tomsk 634050, Russia 

Abstract. Following Ehrenfest's approach, the problem of quantum-classical correspondence can be 
treated in the class of trajectory-coherent functions that approximate as ft —> a quantum-mechanical 
state. This idea leads to a family of systems of ordinary differential equations, called Ehrenfest M-systems 
(AI = 0, 1, 2, . . . ), formally equivalent to the semiclassical approximation for the linear Schrodinger equa- 
tion [16,22]. 

In this paper a similar approach is undertaken for a nonlinear Hartree-type equation with a smooth 
integral kernel. It is demonstrated how quantum characteristics can be retrieved directly from the cor- 
responding Ehrenfest systems, without solving the quantum equation: the semiclassical asymptotics for 
the spectrum are obtained from the rest point solution. One of the key steps is derivation of a modified 
nonlinear superposition principle valid in the class of trajectory-coherent quantum states. 



Introduction 

Semiclassical methods play a distinguished role among asymptotic approaches in linear mathemat- 
ical physics. From the very beginning of quantum mechanics semiclassical approximation has been 
one of the the main technical tools to address its two aspects: pragmatic and philosophical. 

The pragmatic (computational) aspect relies on the presence of a small parameter h as a factor 
next to the derivatives. The pattern is demonstrated in the Schrodinger evolution equation 

ih—^n-^, n = + u(x), p^-ihv^, xew", (0.1) 

ot 2m 
corresponding to the classical Hamilton function 

n{p,x,t) = ^ + U{x). (0.2) 

While Planck's fi. is a dimensional constant, there exists a large class of quantum-mechanical 
problems where a small dimensionless parameter, proportional to h, is present. Accordingly, there 
is a mathematical problem to construct an approximate (with respect to that parameter) solution 
of the quantum mechanical equation. Such an approximate solution is traditionally termed the 
semiclassical asymptotics as h ^ 0. 

The philosophical aspect is related to the correspondence principle, one of the cornerstones of 
quantum mechanics. Despite the fact that quantum mechanics in its axiomatic formalization is a 
self-consistent theory and does not appeal to the classical mechanics, the correspondence principle 
requires the classical equations of motion to emerge from the quantum theory in the limit ft — s- 0. 

Obviously, there is no universal (i.e. physical problem independent) way to obtain arbitrary 
classical values from quantum-mechanical values. In each particular case it is necessary to specify 
in what sense a quantum characteristic becomes classical as ft — > 0. The problem of deriving 
classical equations of motion from those of quantum mechanics in the limit ft ^ is one of the 
principal questions of the quantum-classical correspondence. 

Historically, there are a number of approaches to the problem. One of them is due to Born [1] , in 
which a quantum system is approximately described by the classical statistical ensemble expressed 
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via a semiclassical wave function. A justification of this approach is based on the construction of 
a semiclassical solution to the quantum equation. The time-global version of such a construction 
is known as the Maslov canonical operator [2,3]. In this approach the correspondence principle 
reveals itself in the fact that the principal term of the asymptotic expansion of the quantum density 
matrix is a solution of the classical Liouville equation. 

Another approach, suggested by Ehrenfest [4], is based on the idea that Newtonian equations 
of motion can be obtained in the limit — *■ from equations for mean values of the correspond- 
ing quantum-mechanical observables. More generally, any ordinary differential equations (ODE) 
obtained in the same manner from equations of quantum mechanics can be called classical. The 
correspondence between a quantum observable and its classical counterpart (assuming that such 

exists) is understood as follows: the quantum mean value (^)* of the observable A = A{z, K) cal- 
culated with respect to some special non-stationary states 5'(i; K) must yield in the limit ?i — > the 
corresponding classical observable A evaluated on a certain classical trajectory z{t) in the phase 
space 

Ymi{A)^ ^ A{z{t),Q). (0.3) 

For the Schrodinger equation HO.ll) the Ehrenfest approach relies on states ^{x,t-,h) that are 
localized on the classical trajectory in the following sense: the mean values 

Xk{t,H) = (xk)^ ^ Xk\-^\^dx, pk{t,h) = (pk)^ ^ ^*pk^dx, k^l~fi, (0.4) 

JR" JR" 

of the operators of coordinates x = {xi, . . . , Xn) and momenta p — —ihV calculated with respect 
to such states ^{x, t; h) in the limit h Q 

Xkit) = liinxkit,h), Pfc(t) = limpfe(i,fi), k^T;7i, (0.5) 
obey the classical Hamiltonian system 

mX^P, P=-\7^U{X). (0.6) 

A function ^' for which the limits IjO.SI) exist was called in [5,6] a trajectory- coherent state. 

The technical implementation of the Ehrenfest approach is based on the construction of either 
exact or approximate trajectory-coherent solutions of the Schrodinger equation. Exact trajectory- 
coherent solutions are available only for special Hamiltonians, such as (|0.1I) with quadratic po- 
tential. Examples are well-known coherent and squeezed coherent states [7,8]. An approximate 
[h — > 0) trajectory-coherent solution, called a semiclassically concentrated state can be constructed 
in a much wider class of problems, employing the ideas of the complex WKB-Maslov method [9,10] 
(see also [11-14]). The correspondence principle is manifested in this construction: a trajectory- 
coherent state is an approximate (ft 0) solution of H0.1|) (i.e. it is a semiclassically concentrated 
state) if and only if the trajectory IjO.Sfl satisfies the classical equation 1)0. 6|l . The semiclassically 
concentrated states were first found for particles moving in a potential field [15], and later in an 
arbitrary electromagnetic field [5,6]. Detailed bibliography can be found in the reviews [14, 16]. 

It was found that semiclassically concentrated states exist for linear equations of quantum 
mechanics describing a charged particle with spin or isospin in an external field. In [17-21] the 
semiclassically concentrated states were constructed for the Klein-Gordon and Dirac-Pauli equa- 
tions in an arbitrary electromagnetic field as well as for the Schrodinger and Dirac equations in an 
arbitrary non-abelian field with gauge group SU(2). 

The existence of the semiclassically concentrated states is essential for the approach employed 
in this paper, which consists of the following. Consider an observable A — A(z) whose classical 
analog is A{z). Its mean value in a semiclassically concentrated state can be expressed to any 
accuracy 0(fi(*^+i'/2) via a solution {z{t), A'^{t), A*^(i)} of a finite system of ODEs , 

M 

(i)* = A{z{t)) + Ak{z{t)) ■ A\t) + 0(fi(*'^+i)/2), A'^it) = 0(ft'=/2), (0.7) 

where tensors Ak{z) comprise all partial derivatives of A of order k at the point z, and tensors A*^ 
comprise all moments of order k (see (10 .91) below). 

The dimension of the ODE system is determined by the order of accuracy M. For instance, if 
M = or 1 (it appears that = by construction), we obtain 



{A)^^A{z{t)) + Oih), 
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where z{t) is subject to classical-mechanics equations, in accordance to Ehrenfest's original idea 
HU.3|I . In our approach, by classical equations (of order M ^ 0) corresponding to a quantum 
equa tion we mean that finite system of ODEs whose solution provides accuracy 0(?i'^^^+^^/^) in 
(|U.7|I . and we call it the Ehrenfest M-system. 

The Ehrenfest systems of finite order are truncations of an infinite ODE system, which describes 
evolution of the mean values for a basic infinite set of observables. For the Schrodinger equation 
(|U.1|I the basic set consists of x, p and a special basis of the universal enveloping of the Heisenberg- 
Weyl algebra with generators /, Axk = X k — X k{t), Apk = Pk — Pk, ^ ^ k n, where / is the 
identity operator and Xk , Pk are defined in (|0.4() . The truncations leading to M-systems are made 
due to the estimates for A'^ in (|U.7|) . which allow to disregard within given accuracy 0{h^'^^'^'^'^/'^) all 
variables A*^ for k > M when the means are calculated with resp ect to semiclassically-concentrated 
states. The M — Q truncation is simply Newton's system (|0.t)|) : similarly, the Ehrenfest 0-system 
for the Klein-Gordon equation is the Lorentz equation. 

The infinite ODE system for H0.1|l was derived in [22-24] and it was called the Hamilton- 
Ehrenfest system in [16]. The name reflects a non-trivial fact that the infinite system can be written 
in the Hamiltonian form with respect to a degenerate nonlinear Dirac bracket [25] . A system with 
similar algebraic structure was also derived for the (matrix) Pauli equation [26]. A Hamiltonian 
structure with a degenerate Poisson bracket is also known for the M = 2 truncation [25,27,28]. 
The truncated systems were independently introduced in [29] and used to study quantum problems 
with underlying classically chaotic dynamics. 

In a number of examples this approach was shown to agree with known "classical" equations of 
motion even in the cases where no corresponding classical observables existed. For the Dirac-Pauli 
equation in an external field the Ehrenfest 0-system is a pair of classical equations which are the 
Lorentz equation and the Bargmann-Michel-Telegdi [30] equation in which the field is calculated on 
the trajectories of the Lorentz equation. The order M — 2 truncation obtained in [21] is a Frenkel 
type [31] ODE for spin motion. For the Schrodinger and Dirac equations in external fields with 
gauge group SU{2) the Ehrenfest 2-system [18,20] yields the Wong [32] equation for a non-abelian 
particle with isospin 1/2. More examples of derivations of known "classical" equations from the 
Dirac equation with external fields and the Prock equation are given in [33-35] . 

The Ehrenfest M-system is semiclassically equivalent with accuracy 0(fi,*^*^+^)/^) to the 
Schrodinger equation in the class of trajectory-concentrated states in the sense that it allows us 
to calculate the mean value of an observable directly from the solutions of the system. An explicit 
formula for the state is not required. It was also observed that under certain conditions one 
can obtain asymptotics for pure quantum characteristics, such as energy spectrum series, from 
stationary or periodic solutions of the Ehrenfest 2-system. 

The goal of this paper is to generalize the approach for the case of a nonlinear Hartree-type 
equation. In particular, we consider the following equation of self-consistent field 

at 

n.m = -^Vl + U{x) + ^ jv{x,y)\^>{y,t)\^dy, feR", ^^'^^ 

where U{x) and V{x, y) are given smooth potentials of the external electromagnetic field and the 
self-consistent field respectively, and k \s a. constant. 

There are at least two reasons why the problem of quantum-classical correspondence was not 
considered in the spirit of Ehrenfest's approach, neither for a nonlinear self-consistent field (|0.8|l . 
nor for more general Hartree-type equations. 

First, the operator 7ijj(^') does not have a natural classical analog in the traditional sense, thus 
it is not a priori obvious which dynamical ^'-independent system is an appropriate candidate for 
the "classical" system in the limit Q. 

Second, it is not clear whether the nonlinear quantum equation has either exact or approximate 
(as ?i — *■ 0) solutions that are trajectory-coherent in the sense of H0.5|l . 

In the framework of our approach, a solution of the correspondence problem includes three 
stages. 

1. First, for an arbitrary one parameter family of phase space trajectories Z{t, h) = {P{t, h),X{t, h)), 
t € IS., we introduce a class V\{Z{t,h)) of trajectory- coherent functions. Exact construction is 
given in Sect. 2. 

Let a, /3 e Z" be multi- indices, |q;| = T%=i '^k, and X" — YYj=i -^'j' ■ Let A^^ be an operator 
with Weyl symbol Aq^ (p, a?) = {x — X{t))°' {p— P{t))^ . The following property is essential for 
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further steps: the centered moments 



Aa^(t,?i) = ^A„^^^, (0.9) 
calculated with respect to functions from class 'P^{Z{t, h)) satisfy the estimate 

Consequently, k-th order moments (i.e. those with \a\ + = k) are Oih^/'^). 

Next, we assume that equation (|0.8|) has either an exact or approximate (with accuracy 

0(/i^*^+^)/^), M ^ 0) solution ^' in the class of trajectory-coherent functions. Using an ap- 
proach similar to the linear case, we derive an infinite Ehrenfest system for the nonlinear 

equation and its finite M-truncations for {{P{t, K), X{t, K)) , A„/3(i, h), \a\ + |/3| < M}. 
Details are provided in Sect. 3. In particular, the principal {M = 0) Ehrenfest system for the 
nonlinear equation of self-consistent field has the form 

mX = P, 

(0.10) 



y = S=X 



Note that when xr = 0, this system turns to (|0.6() . Similarly to the latter, the system (|0.10|l 
describes mod h}^"^ the trajectory where the trajectory-coherent solution €V^is localized. 

3. Given a semiclassically-concentrated solution of equation (|0.8|) with accuracy 0{h'^^^^^^^'^), one 
can obtain a corresponding solution to the Ehrenfest system of order M by evaluating mean 

values of operators x, p, A^^ with respect to that solution. Our goal is to show that it works 

the other way around as well: quantum characteristics can be found with accuracy 0(fi,*^*^+^)/^) 
from a solution of the Ehrenfest i\f-system. 

In Section 4 we study the Ehrenfest 2-system for the Hartree-type equation and in particular 
its solutions corresponding to the rest point of the classical (M = 0) system. Based on this 
calculation, in Section 5 we reconstruct the asymptotics of the energy spectrum for the Hartree- 
type equation. Examples in Section 6 illustrate general results. 

The key point of the whole approach is to obtain quantum characteristics without explicitly 
solving the quantum mechanical equation. An explicit formula for the solutions is not required 
anywhere in the derivation. A basic assumption of the paper is the existence of a semiclassically 
concentrated solution of the Hartree-type equation in a class V^. This assumption can be justified 
by explicit construction of formal asymptotic solutions using finite dimensional Ehrenfest systems 
[36-40]. 



1. Hartree-type equation 

By the Hartree-type equation we mean the following equation 

{-ihdt+'HU'^)}'if = 0, 'H^{^)^H + hV{-^), *eL2(M;'). (1.1) 

Here 

n - H(i), V{^) = J dy^*{y, t)Viz, w)^{y, t), (1.2) 

where the pseudo-differential operators 7i(z) and V{z,w) with symbols H{z) and V{z,w) respec- 
tively are functions of non-commutative operators 

z=(-iWa;,x), w ^ {-ihX7y,y), x, j/eM", 

Function ^P* is the complex conjugate to 5*, >«: is a real parameter, ft > is a small parameter. 
The operators z and w satisfy the following commutation relations 



[zk,Zj] = [wk,Wj] = ihJkj, 
[zk,Wj]=0, k,j = l,2n, 
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(1.3) 
(1.4) 



where J — \\Jkj\\2nx2n is the standard symplectic matrix 




and [A, 13] — AB — 13 A denotes the commutator of A and B. 

In this paper all functions of non-commutative operators are Weyl-ordered [41,42]. The action 
of the operator Ti in this case can be written as 

■H^{x,t,h)^j^^^ J dydpexp(^^{x-y,p)^n{p,^^)^{y,t,h), (1.5) 

where Ti.{z) — TL{p, x) is the Weyl symbol of the operator TL, and {x,p) = X]fc=i ^kPk- 

Remark. In the particular case when the Weyl symbols of operators 7i(z) and V{z,w) in H1.2f) 
have the form 

^2 

2m 

the equation (|l.l|l yields the equation of the self-consistent field in the form (|0.8|l . This differential 
equation with integral nonlinearity plays a fundamental role in quantum theory and nonlinear 
optics [45,46] and in the theory of the Bose- Einstein condensate [47]. In the latter the solution 5* 
represents the wave function of the condensate, while the non-local potential V{x, y) describes the 
interaction of condensate's particles with external field. 

In this paper we deal with asymptotic solutions of the equation localized in the usual 

mathematical sense rather than in the sense of (|0.t)|) : namely, functions or formal series '^{x,t,h) 
must belong to the Schwartz space with respect to the variables x G M". We require the Weyl 
symbols Hiz) and V{z, w) of the operators H and V{z, w) in H1.2|l to belong to one of the 
classes [3, p. 13]: they must be smooth functions of at most polynomial growth with all derivatives, 
such that the following conditions hold. 

Assumption 1. The functions 'H{z) and V{z,w) are infinitely differentiahle for all z G M^" and 

w £ R^", and for any multi-indices a, /i G Z^" there exist constants Ca, Cap, and m ^ such that 

i^Cap{i + \z\r{i + \w\r. 

The notations here are as follows: 

ai + a2 + ■ ■ ■ + a2n, 

(1.6) 

dz'^^dzr...dz^i-- 

Note that for our method it is essential to have smooth symbols Ti-iz) and V{z,'w). Asymptotics 
for Hartee-type equations with singularities are a subject of a number of publications, see e.g. [43] 
and references therein. 

Now we introduce a vector space in which asymptotic solutions to the equation Hl.lfl will be 
sought. 



a<^'n(z) 



^C„(l-t-lzl)^ 



Oi" ' V [Z,W) 



a = {ai,a2, ■ 



■ , a2n), a.j ^ 0, \a\ = 



2. Class of trajectory-coherent functions 



We will construct asymptotic solutions of equation Hl.l|) with the following features: they have a 
form of generalized solitary waves and singularly depend of the small parameter ft ^ 0. 

Such a solution relies on a phase space trajectory Z{t, h) = {P{t, h), X{t, h)) and is trajectory- 
coherent in the sense of (|(J.5|) . (|0.4|l . We denote the class of trajectory-concentrated functions by 
V^{Z{t, h)) and define it more precisely as 



vi^rl{z{t,h)) = \<i>:'fix,t,h) 



exp 



-{S{t,h) + (Pit,h),Ax)) 



(2.1) 



where function Lp{^,t,h) belongs to the Schwartz space § with respect to variables ^ G M", is a 
smooth function of t, and regularly depends on ft as ft —> (the term function is used throughout 
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in the sequel, although ip{(^,t,h) may be in fact a formal series in powers of h^^^). Here Ax — 

X{t,h). The real function S{t,h) and 2n- vector- f unct ion Z{t,h) also regularly depend on h 
When an asymptotic solution of equation is being constructed, these functions, 



0. 



X 

as h 

as well as the amplitude t, h) are to be determined. Short notation Pf^ for Vf^{Z{t, h)) will be 
used when it does not lead to confusion, and || • || will denote the L2-iiorm for functions from V^. 
In the expression for the norm of a function from class 7^^ the argument t may be omitted, and we 
write 

It will be shown in Section 3 that the functions h) and S{t, K) are uniquely determined 
by the Ehrenfest system corresponding to the Hamiltonian of equation In the linear case 

{x = 0) the vector-function Z{t, 0) and the scalar function S{t, 0), defined by Hamiltonian function 
Ti.{p, x), are the classical-mechanics phase space trajectory and the classical action respectively. As 
an example of amplitude, dynamical coherent states for quadratic Hamiltonians in the form of 
Gaussian function can be given: 

V.(f,t) = exp[^(e,Q(t)O]/(0, 

where Q(t) is a complex symmetric matrix with positive imaginary part, and the time dependent 
factor f{t) is given by 

t 

fit) = VlmQ(t)exp[-^y"lm g(r)dT 



(see for details [16]). 

Consider important properties of functions from class 7^^. Their proofs are given in [16]. We 
briefly reproduce some of them in Appendix B. 

1. Let $ belong to the class 'Pf^{Z{t, h)). Introduce operator {Az}" with Weyl symbol (Az)" = 
(Azi)"i •...•(Az2„)"^", and 

Az = z - Z{t,h) = {Ap,Ax), Ap = p- P{t,h), Ax = X - X{t,H). 

ations for mc 



Then the following asymptotic estimations for moments Aa{t, h) of order |a|, a e Z^" hold 



(2.2) 



Denote by 0{K^) an operator F such that for any function $, from the class Vl^{Z{t, h)), the 
asymptotic estimate holds 

IIF$|| 

^ = o(r), ;.^o. 

2. The following asymptotic formula holds 

{Az}" = d(nl"l/2), aeZf, h^O, (2.3) 

in particular 

{Axfe} = 0(fti/2), {Ap,} = d{h'/^), k,j=T;7i. (2.4) 

3. For functions $(f, t,ft) G V^{Z{t,hj) the following limits hold 

hm ^|$(f , t, h)f = 5(f - X{t, 0)), (2.5) 
}^^,7;^\mt,fi)f ^ 5{p- P{t,0)), (2.6) 

h^o ||$|j2 

where $(p, t, h) = Ffi^^^j;^{x, t, H), Fj^^^^j; is ft^^-Fourier transform [3]. 

Denote by {A(t))^ the mean value of a self-adjoint in L2(K") operator A(t), t e R^, calculated 
with respect to the function ^(a;, t, h) G V^. Then we have 
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4. For a function ^{x, t, h) G V\{Z{t^ h)) and an operator A{t, K) with Weyl symbol A{z, t, h) 
satisfying the first inequality in Assumption 1, the following equality holds 



]ira{A{t,h))^ = lim^j^{<l>{x,t,n)\A{t,h)\^x,t,h)) ^ A{Z{t,0),t,0). (2.7) 

The limiting nature of conditions H2.5(l . (|2.6(l and asymptotic character of estimations (I2.2|l - 
H2.4|l . holding in the class of trajectory-concentrated functions allows the construction of an 
approximate solution ^Pas — "^stsix, t,h) of the Hartree-type equation for any finite time interval 
[0, T], in the following sense 



-in^, + n + >€V{^, 

at 



^-as = 0{K>), (2.8) 

^^,&vUz{t,h),S{t,h)), ie[o,r], (2.9) 



where 0{H'^) denotes a function g^'^'>{x,t, K) which represents the error for equation and the 

error obeys the estimate 

^rnax^ \\g'^'^\x,t, h)\\ = 0(^9), ft ^ 0. (2.10) 

Following paper [16] and having in mind properties (|2.5|) . (|2.6() . we call such a functi on ^ as(a?, t, h) 
a semiclassically-concentrated solution (modft", H —> 0) for a Hartree-type equation (|l.l(l . 

The semiclassically concentrated (modft") solution ^'(^^(x, ft) of the Hartree-type equation 
is a formal asymptotic solution evolving from the initial state {x, ft) chosen in the class of 
trajectory-concentrated functions V^(zo, Sq). Here zq = (j)o,xo) is an arbitrary point of the phase 
space Mp", and the constant So = S{0,h) can be set equal to zero without loss of generahty We 
will denote the class of initial {t — 0) trajectory-concentrated functions by V^{zo). Exphcitly, 

po(zo) = |v:V'(x,ft) = ^o(^^,fi)exp{^(po,f-xo)}, ^o(f, ft) € §(R^)|, (2.11) 

Let us give two important examples of the amplitude function in the initial state (|2.11fl . 
First, 

^o(0=e-<^•^^>/^ 

where the n x n-matrix A is real symmetric and positive definite. In this case relation (|2.11|) 
describes the Gaussian wave packet. 
Second, 

where the n x n-matrix Q is complex symmetric and has positive definite imaginary part Im Q, and 
ly — {vi, . . . , i^n) is a multi- index of the multi-dimensional Hermite polynomial Hi,{ff), ff e M" [48]. 
In this case ip G T^hi^o) defines the Fock state of a multi- dimensional oscillator. 

A construction of a semiclassically-concentrated solution (modft"), of the problem (|l.l|l with 
initial state from (|2.11|) is based on the solution of the Ehrenfest system, to which we turn our 
attention now. 



3. The Ehrenfest system of equations 

Let symbols Ti.{z) V{z^w) satisfy Assumption 1. Then operator H(z] 11.211 is self-adjoint with 
respect to inner product (^['l') in space L2{K^), and operator V{z,w) (|1.2() is self-adjoint for the 
inner product in space L2(K^y). Thus the norm of the exact solutions of (|l.l|l is preserved by time 

evolution: ||^'(t)|| = ||*o||- The mean value (A) = {A)q, = (^lAI^") of an operator A{t) = A{z,t), 
calculated with respect to these solutions obeys 

l(A{t)) = + ^{[n{z),Ait)]) + '-^{1 dy^*iy,t^h)[Vi^^ (3.1) 

as an implication of the Heisenberg equation for evolution of operators. Equation 1)3. l|l is called 

the Ehrenfest equation for operator A(t) and function ^'(af, t, ft). Our choice of this terminology is 
justified by analogy with the linear case (>£■ = 0) in which equation Hl.l|) becomes the Schrodinger 
equation, while equation H3.1|l is the Ehrenfest equation [4]. 
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To derive the Ehrenfest system from the Ehrenfest equation (|3.1() we take for 4" a solution of the 
Hartree-type equation in the class of trajectory-concentrated functions, and for A operators 
z = {p, x) and {Az}" with Weyl symbols {Az}", a <E Z^", where 



Az = z- z{t,n), Z{t,h) = {■^{t)\z\^{t)). 
We represent operators 7i(i,t) and V{z,w,t) in the series form 



OC 

n{z) = n{z{t,h))+ E -nf,{z{t,h)){^^zr, 

\^,\=l 

oo oo . 

viz,w)^vizit,h),zit,h))+ J2 E —v,Azit,h))A,{Azr, 



(3.2) 



(3.3) 



\,y\ = lM=2 ^ 

d\^\n{z,t) d\^^+-\v{z,w,t) 

^A'l^'*)^ ' V^u{z,t)- 



dzt^dw'^ 



where A^ as defined in l|2.2(l . Thus to derive the system we need to evaluate commutators 
[zk,{Az}i^] and [{Az}", {Az}^] for k = and 1^*1 ^ 1, \v\ ^ 1. This has been done in the 

linear case {k = 0) [16,22-24]) using the formula for composition of Weyl symbols A{z) and B{z) 

(see e.g. Appendix in [41]), defining the symbol C(z) of the product of operators C — AB 



2 ih d 



(3.4) 



Here the number over an operator refers to the order of its action onto the target function. 

This way we obtain an infinite system of equations for Z{t,h), Aa{t,h). Keeping only the 
moments up to order N, we obtain the following finite system of equations 



^1 ^1 
E + E ^V,^,{z)A^A,y 



1/^1=0 



N 



E 

1^+71=0 



_ . x l7l-i [(-1)1^^1 - {-l)\f-\]a\/3W{a - -f)e{(3 - 7) 
/ 7!(a-7)!(/3-7)!Ai! 



(3.5) 



N 

(H^(z) + i< E -,y^.u{z)A,)A 



2n 

a-7 + J,3-J7 - E Zk<^kAa(k) 
k=l 



with initial conditions 

z|^^Q = zo = (VI^IV'), A^I^^Q = (V-liz-zo}"!^), aeZ2«, |a| ^ iV. 
Here >c = K\\ip{x,h)\\'^ , and 'ip{x, H) is the initial state from 7')?'(zo) H2.11|l . 

2n 



(3.6) 



d\^\nM ,r ai^+Hi-,(z,zz;) 
n,f,{z) = — — — , v,f,^{z) = 



dzt^dw'^ 



eia~(3) = Y[9{ak-Pk), (3.7) 



k=l 



a^{ap,ax), Ja^{ax,ap), a{k) ^ {ai ~ 6ki, ■ ■ ■ ,a2n - 5k2n)- 

As in the linear case (x = 0) (see [16]) the system (|3.5|l will be called Ehrenfest system of order 
N. Due to estimates (|2.2ll this system is equivalent in the class of trajectory-concentrated states 
to the Hartree-type equation 11. 1|) with precision 0(fi,^^+^^/^). 

Introduce notations 



(3.8) 



2nx2n 



i5(z, w) 


^n{z) + kv{z,w), 


^z{z) = 


^z{z 






W — Z 


(9^io(z, w) 
dzjdzk 


W — Z 





d9){z, w) 






10 — z 


dzj 


W — Z 


•1 

lx2n 



Then for = the Ehrenfest system (|3.5|l has the form 

z = JSjz{z), 



(3.9) 



and for = 2 we obtain 



A2 = Jjo..(z)A2-A2io,.(z)J, AT=A2, 
where A2 is a symmetric (2?! x 2n)-matrix of the second moments (A2 = |1 Aji||2nx2n) 



(3.10) 



A2(0 



[t) (Tpa;{t) 
(t) (Jxx{t) 



(3.11) 



\(ApkApi))\\ 



\{AxkAxi) 



-{AxkApi + ApiAxk] 



with (71 X n)-blocks 

(^xp{t) — \WxkPi{^)\\nxi 

System can be written in an equivalent form introducing matrix A{t) via relation 

A2{t)^A{t)A2{0)A+{t). 

Then the system becomes 

1 + AA^A+a,) + ]^{d^,AAlA+d^)\^^{z,w] 




A(0) = L 



(3.12) 



Note that the initial state function can be excluded from the initial conditions (|3.(j|) for the 
system H3.1()ll if they satisfy an infinite system of ineq ualities (generalized Heisenberg uncertainty 
conditions) and for t — sat isfy estimations (I2.2|l . Obviously, all Heisenberg inequalities are 
consistent with equations (|3.5(l . 

The uncertainty relation corresponding to the Ehrenfest system (|3.1Uf) can be rewritten as the 

ih 

condition that matrix A2{t) + — J is positive definite [50] (see also [16,51]). 

The fact that equation (|3.5|) is a finite system of ordinary differential equations for functions 
z, Aa equivalent to problem p.l|l with precision 0( ^*^^+ -'-^/^), suggests that there is a mechanical 
system with finite degrees of freedom described by H3.5|l . and thus an approximate semiclassical 
model of a Hartree-type equation allows an exact interpretation in frame of classical mechanics. 
The number of degrees of freedom of such a mechanical system grows with precision N . A study 
of such classical systems with quantum origin by methods of classical mechanics constitutes a 
separate interesting direction of research [27,28]. 

4. The Ehrenfest system: 

the small dispersions approximation 



Denote by q operators 



A^' = -{AzkAzi + AziAzu); j,k,l = l,2n,fc ^ / 



(4.1) 



Their mean values in state ^'(x, t) with initial condition ^'(x, t) = ip{x), are denoted by 



Q^{t, h) - {-^{x, m\^{x, t)) = z{t, h, sl),A2{t, h, . 



(4.2) 



The mean values g^{t,h) obey Ehrenfest system (|3.10() with initial condition g^(0,fi.) = 0° = 
(■0101-0). When it does not lead to confusion the explicit dependence of the initial fimction will be 
omitted in the notations. 



*For |a + /3[ ^4 the system of inequalities is well known (see e.g. [51], and the bibliography there). 
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Estimations (|2.2|l suggest that the Ehrenfest system can be solved approximatly with respect 
to a small parameter H using the expansions 

g{t, h) = (t, h) + ?i0(i) {t,h) + ..., (4.3) 

or 

Z(t) = Z(0)(i) + ?iZ(i)(t) + ..., 

(4.4) 

/^2{t) = /^f\t,h) + h^'^^\t,h) + .... 

Here Q^'^\t,h) = [z^^\t), lS.f\t,K)^ is the principal term of the solution of system l|TTn|) . Sub- 
stituting (|4.3|l into system l|3.10ll we obtain equations for z*^"' = Z^^\ z'^' ~ Z'^^\ and with 
precision 0{tv'^^) 

^ j^4z(")), 

i^'^ = J^,M°^V^^ + Fiz(°\A^°^), (4.5) 

A^°)= Ji0,.(z(o))Af -Afi3..(z(°))J. 
Here 

F{z,A2) = ^Jd,Sp\\sj,,{z,w) + KV^^{z,w)]A2} . (4.6) 

The first equation of system (|4.5|l coincides with H3.9|l and is similar to the classical-mechanical 
Hamilton system in the linear case, however, in the nonlinear case (i< 7^ 0) the system is not 
Hamiltonian. 

Consider the following auxiliary system of equations which can be referred as the pseudo- 
system-in-variations for the solution Z^^\t) 

dfc = Ji5,,(z(")(i))afe k^T;7l, (4.7) 

with normalization condition 

{ak{t),ai{t)} = {al{t),a:{t)} = 0, {aUt), ai{t)} ^ ~2iSki, (4.8) 
where {w, u} is a skew symmetric inner product in M^" 

{v, u} - (v, J^u) = Wl n - 1^6: (4.9) 



Ya ) V n 

If a solution z*^") = Z(°^(t) as well as a complete set of solutions a^it) of the pseudo-system-in- 
variation are known then the general solution of the two last equations in H4.5|) has the form 



z^'-\t) = E + . (4.10) 
fe=i 

Af^{t)= A{t)VA'{t), (4.11) 
where scalar functions hk(t) and the 2n x 2n-matrix A(i) are as follows 



t 

1 



&fcW = -^ J{aUt),Fit)}dt + b,, F(t)=^^(Z(")(t),Af (i)), (4.12) 


A(t) = (^a,{t),a2{t),...,an{t),al{t), ■■■,<{*))■ (4-13) 

Here 6^ are constants of integration, V is an arbitrary 2nx2n constant matrix and F{Z'^'^'> (t), A'^^ (t)) 
is defined in H4.6|l . Thus in this approximation the total solution is determined by solutions of the 
modified classical system (|3.9ll and pseudo-system-in- variations l|4.7|l . 
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5. Energy spectrum and the Ehrenfest system for quantum 
means. 

Consider the problem of defining the energy spectrum for a Hartree-type Hamiltonian Ti^^ from 
the dynamical Ehrenfest system for quantum means. Denote by S[gj^(t, h)] the generalized action 
[36,37] along the trajectory g-,j,{t,h) H4.2|l of the Ehrenfest system H3.1U|I 



S 



dt 



P(i,0^),X(i,0°))-io(z,u;) -|Sp(K,^(z,u.)A2(t,4) 



-ft(p(")(t,0°),xW(t,0°)) + %o,X^'^) 
Substituting (|4.4|l in H5.1|l we get 



w=z=Z{t.,sl) 

(5.1) 



S 



q^M] =y'rfi[(p(o)(i,0O),i(o)(<,gO))-f,(z,™)- 



-hk{V^{z, w), [t, 4)) - |Sp(k,»(z, HA^") {t, gO ) 



^=2=Z(")(t,s") 



+ C»(ft3/2). (5.2) 



The connection between the energy spectrum and the time-dependent solution (t, K) is es- 
tablished as follows 



Statement 5.1. Let the Hartree-type stationary equation 



(5.3) 



have pure discrete non- degenerate spectrum, and functions Lp^, form a basis in T'^{zq) (|2.11() . Let 
a solution g(t, h) of system ()3.1U|) exist in the class of quasi-periodic functions [49]: 



.{h)e^' 



(5.4) 



and satisfy the generalized Heisenberg inequalities. Here Q^^' (ft) are generalized Fourier coefficients 
for function 2{t,h). Then the energy levels and the frequencies oj^v' are connected by the 
following relation (mod ft^/^) 



hiu,,, ^E,- E,, + -^S[Q^^{t,h)] - i5[g^,^,(i,ft)], uu' e II] 



(5.5) 



Proof. At the initial moment of time function g(t, H) satisfies the generalized Heisenberg 
inequalities. Thus there exists a function ip € 'P^{zq) such that with the accuracy 0{h^/^), g(0, K) = 
0^(0, ft) = (see [16]). 

Since vectors ipi, form a complete set and belong to 7^)^(2:0)1 every semiclassically-concentrated 
solution of p.l|l with the same Hamiltonian as in (|5.3fl . and initial state -0 G T^ni^o): can be 
represented in the form (see Appendix A): 



oc 

^(x,i) =exp(^5[0^(t,ft)]) Y aexp(--^S,t--^5[g^„(t,ft)])(p.(f) + 0(ft3/2). 



(5.6) 



Using function '^{x^t) in this form for evaluation of the mean value of operators g H4.1|l and 
comparing with (|5.4(l we get 1)5. 5|l , which completes the proof. 

If the ground level of energy Eq as well as the complete set of frequencies LO,jyi are given, then 
the problem of reconstruction of the entire spectrum Ei, is a well known problem in spectroscopy. 

The invariant manifolds of classical Hamiltonian systems are known to be important for clas- 
sification of spectral series in the linear theory (>r = 0) [10]. 

We are able to solve the problem in the semiclassical approximation using invariant manifolds 
g(°)(i) of system (|4.5|l for classification of the spectral series. 
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In the simpl est c ase when the invariant set consists of a single rest point z'^^-' = = const 

it follows from 14.511 that 

f),{z^°^)^dAn{z)^kV{z,w)\ =0. (5.7) 

Identify the rest point Z'"' with point zq which defines class Vf^^zo) H2.11|l . In other words we take 

initial states localized at the rest point zq = Z^"^. 
We will proceed under the following assumption. 

Assumption 2. Let the symplectic (2n x 2n)-matrix JS)zz{zo) evaluated at the rest point zq have 
n distinct pure imaginary eigenvalues iflk, > 0, fc = l,n (and n eigenvalues complex conjugate 
to them, ~i^k, k = l,n). 

In the linear theory, Assumption 2 implies stability of the rest point in the linear approximation 

[9]. _ 

Note that under Assumption 2, solutions of the pseudo-system- in- variation H4.7|l have the form 

afc(t) = exp(irifct)/fc, k=l,n, (5.8) 
where fk is the eigenvector of the pseudo-matrix-in-variations, evaluated at the rest point 

J^zz{zo)fk = i^^kfk, ^k ^ ^j, j, k=l,n. (5.9) 

The eigenvectors fk,k = l,n are normalized, without loss of generality, by condition H4.8|l . 

Statement 5.2. Under Assumption 2, the energy spectrum of the Hartree-type operator (|5.3|l can 
be found as 

E,^S)(zo,zo) + hJ2hk(^iyk + ^^ +0{h^/^), (5.10) 
where multi-index u — (i^i, . . . , z/^^) has all non-negative components, 9}{z^ w) is defined in H3.8fl and 
flk^flk + -^(/fc,K,«,(zo,zo)/fc) + 

n 

+ ReY,i^{yUzo,zo)Jk){f^,d,){f*,[^,,{z,w) + kV^Uz,w)]fj} _ _ . (5.11) 



2n, 



Z—W — Zq 



For 1/ = our result 15.1011 coincides with results obtained in [52] with precision 0(fi.'^/^). 
To prove Statement 15 . 21 we need several Lemmas. 

Lemma 5.1. Under Assumption 2, the spectrum of Hartree-type operator (|5.3|l is equidistant 
(modh^^^) and can be found as 

n 

E, = Eo + hY,^ki^k, (5.12) 

fe=i 

where Eq =const is the ground energy level and flk are defined in (|5.11(l . 

Proof. As the first s tep, we establis h the linea r rela tion between fTfe, k ~ 1, n and lui,^' in the 
Statement 5.1. Using formulas and with (Ol) . (j^^ become 

n 

Z(^\t) = ^Rc [{f^,Sk{t)}fk + fefce'^^Vfe] , bk = const, 

''=n^ (5.13) 

^2°^W= E 2(^^"^'^ + ^^^')^^''^^'''"'''^'' = I?,/=const. 
j.i=i 

Here the 2n- vector ^k{t) is defined by the relation 

dkit)^^Jd,Sp{\^),,iz,w) + ?cV^Uz,w)\j^k{t)} , (5.14) 

Zfl I L J J w—z—ZQ 
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where the (2n x 2n)-matrix Tkit) has the foUowing structure 

n 

Using the estimate (I2.2f) . solution of the Ehrenfest system (|3.5|) can be written in the form 

Q{t,h) = (^.Q + hz^^\t),i^f\t)^ +o{h^/^). (5.16) 

Introduce vectors 

/2== (/ii,...,^„), = {i^i, . . . n ^ {ni, . . . ,nn), (5.17) 

and le t denote the sum of the absolute values of vector's components. Substitute H5.13|l into 
(|5.16|l and rewrite it as 

0(t, h)^J2 B^me'^^'"^^' + 0{K"^). (5.18) 

The right hand side in H5.18|l is an almost periodic function with n frequencies f2i,...51„. The 
exponents in H5.18|l are linear combinations of the frequencies fc = 1, 71 with integer coefh cient s 
/Life, and thus linearly depend on the multi-index fi € Z". Consequently, frequencies uJw' (|5.4|I . 
must be linear functions of v and v' . Taking into account that by interchanging v and v' frequencies 
ijJvv' H5.5|l change sign, and by preceding arguments, we get 

uj,,' = {u-u'M) ^{t^A. Ia?K2. (5.19) 



The step 2 of our proof consists in recalculation of the right hand side of the formula H5.5|l . 
Energy level E^, with precision 0{hP^^) can be found from the mean value of Tij^ in the stationary 
state (fii, e V^{zo) as follows 



1 



(5.20) 



(5.21) 



+ - Sp { [io,,(zo) + i<K,^(zo, zo)] A^} + 0{h^/^). 

Here the equality (|5.7|) was used. 

Note that for the stationary state (pi, we have for coefficients in (I5.13f) 6^ = 
since the functions are time-independent 

k=l 1=1 

Here vector is defined by 15.14|l with 

n 

Substituting formulas (j5.21|l into H5.2QI) we get after some calculations using ((5 

n 

E, = io(zo, zo) + MK«(zo, zo), Zj"') + ^ [rife + |(/*, y^„(zo, zo)/fc)]l?["^ + 0{fi"^){^.2?,) 

fe=i 

Similarly, from H5.2|l we find 

1 " 

T^[0^J = -^(^0,^0) - MK,(^o,^o),^l'^) - ^<5I<-^fe'^»»(^0'^o)/'')^i'^^' (5-24) 



(5.22) 



fc=i 
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Now substitute (jF^^ . lfOljl in to find that 



From this equation along with (|5.19|) we can find D^'^'' in terms of uf'' as follows 



V 



(5.25) 



Constants v'^p will be defined later. 



To illustrate how the equation H5.25|l was obtained we temporarily consider case n — 2, without 
loss of generality. Then condition ^ 2 in H5.19|l reduces our consideration effectively to the 
three cases. 

Case 1. If j-ii = 1, /i2 = then Lijy,ji ~ f^i and ui — v'l = 1, 1^2 ^ 1^2 = 0. Consequently, 



hQ.1 = Q.i{D\ 



) + n2iD. 



(l'l,l'2) 



D. 



and thus 

Case 2. If /^i = 0, /i2 = 1 then cj^j/' = il2 and j^i — i/J = 0, 1^2 — = 1- Consequently, 



(''1, 1^2) 



(1^1,1/2-1)^ 



n2(L'; 



(''1, 1^2) 



(1^1,1/2-1) 



and thus 

'1 — -^1 1 ^2 ~ ^2 

Case 3. If //i = 1, /i2 = 1 then cj^i^' = fii + ri2 and i^i — = 1, 1^2 ^ ^^2 = 1- Consequently, 



{1^1,^2) 



and thus 



(1/1-1,1/2-1) 



h, D 



{vi,V2) 



= D 



(l/l:l'2) 

2 

(l/l-l, 2/2-1) 



^(.1-1,.2-1))^ 



Note that the third case result can be decomposed into the first t wo. It is apparent that similar 
derivations are possible in any dimension n ^ 2, and that equation (|5.25|l holds. 

As the third step of the proof, we finally obtain the statement of Lemma 5.1. Using H5.21|l for 
Zi^^ with if?:^ for 5-^, from g^J and ifrT^ . we find 



(K,(^o, z^),Z^l'^) - Re ^ 5] ^(^-(^' h 



k=ij=i 3 
From H5.23|l . using and (|??^ we get 



n 



(5.26) 



(5.27) 



where fij is given by H5.11|l . Finally, denote by 



(0) 



(5.28) 



to obtain (|5.12|) . Thus the lemma is proved. 

Lemma 5.2. Under Assumption 2, the ground energy level for the Hartree-type operator H5.3|l is 
given (mod h'^/'^) by 



Eq = 9]{zq, ^o) + -^^ ^ ^j- 



(5.29) 
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Proof. The ground energy level Eq is given by (|5.28() . which follows from (|5.27() for v — 0. Con- 
stants 2?^*^^ will be chosen in order to minimize the uncertainty condition and obey the quantization 
condition (|5.5|) . This idea is supported by physical models including the harmonic oscillator and 
the Coulomb potential [53,54] (see [51] and references there). 

ih . 

As we have mentioned above, matrix A2(t) + — J is positive definite. Thus for any vector v 
we have 

V+ A2(t) + yJ f ^ 0, 

where equality corresponds to minimization of the uncertainty relation. 

Take v — Jfj,j = l,n. Then using 15.21|l with v = and orthogonal condition (|4.8I) (which is 
valid for vectors fk, as well as for ak{t)) we have 



(0) 



So, choosing 



(0) _ 



(5.30) 



we minimize the uncertainty relation, and obt ain 15.2911 . 

Then formula H5.10|l directly follows from (|5.12l) and H5.29|l . Thus statement 15. 21 is proved. 



6. Spectrum for the oscillator with nonlinear Gaussian po- 
tential 

6.1. Spectrum for the oscillator in constant magnetic field and nonlinear 
Gaussian potential. 

In this section we illustrate the method described above with an example of a Hartree-type equation 
p.l|l whose linear part corresponds to an oscillator in a constant magnetic field, while the nonlinear 

part is described via Gaussian potential. The linear part Ti. has the form 

n = -Li0-lA{s)r + '^x\ (6.1) 

2m c 2 

The external field in the operator H6.1|l is a superposition of a constant magnetic field H = (0, 0, H) 

^ 1 ^ k 
with vector potential A = —Hxx and an oscillatory field with scalar potential ' "^^"^ non-local 

operator V{'i') in (|1.2() has the form 

{x-yf 



V{^)= V{x,y)\^fiy,t)\'dy, V{x,y)^Voe^p 



272 



(6.2) 



Here H, Vq, fc, 7, e, c are real parameters of the model. 

We will be using notations ujh, ujq, and uj^i for the cyclotron frequency, oscillator frequency, 
and nonlinear frequency respectively 



eH k \kVo\ 

^^H = , Wo = -\/— , CJni = W -, (6.3) 

mc \ m y m7"^ 



where >c — xll'i'W . We also introduce 



-0^/1+ (|^]'. (6.4) 



To construct a solution 5" e Vl{Z{t,h)) of equation Hl.l|l . H6.1|l - H6.2|l . we are using a phase 
space trajectory Z{t,h) = {P{t,h),X{t,h)) which obeys H0.5|l . (|0.4() . 



15 



Using notations H3.8() we have 



Hxi + x^) 



2 

/ 1 



Pi 



m 
1 



-Pi 



X2 



^Xi + ^iP^X.-PM; 
\ 



(6.5) 



2 
2 

m 



^P2 + muj^Xi 

"^P.+mujlX, 



\ mujQXz 



(6.6) 



where z = (P,X), P = {Pi,P2,Pz), X = {Xi,X2,Xz), and notations (jHS, are used. The 

matrix of the second derivatives becomes 



^zz{z) 

^pp{z) 
^xx{z) 



■^pp(-^) ^px (^) 
S^xp{z) Sjxxiz) 



1 1 1 

2 



I^PfcPi(2)ll3x3 = diag[ 

Il^)x,a;,(z)||3x3 = diag(m(w^ - Tjujli),m.{ujl - 77w^i):"^(^o - W^i)); 



/ 



%:r(^) = ||^)pfc:r,(2:)||3x3 





2 

\ / 



(6.7) 

(6.8) 
(6.9) 

(6.10) 



Here r\ — sign(iirVb). 

Recall, that to find the spectrum corresponding to the Hamiltonian in we need only 

bounded solutions of the Ehrenfest system. The first equation of system (|4.5|l describes Zo{t,h), 
and can be integrated independently from the other equations of the system. The last equation of 
14.5|l describes the second moments A2(t) and depends on the solution of the first one. Therefore 
we start with solving the first equation of H4.5|l . The simplest stationary solution of it is the zero 
solution 

Zo{h) = (Po(ft), ^o(^)) ^ = (0, 0, 0, 0, 0, Oy. (6.11) 

Then the corresponding eigenvalue problem (|5.9() has solutions ili = u!+, fli = co-, f^a = ^s, 
where the Ritz frequencies 



UJ+ 



and the eigenvectors are 



2 ' 



/i 

h 

h 



9 2 '^a 



2 ' 



1 i 1 

^(5o,«5o,0, , — ,0)T, 

V2 50 go 

1 i I 

— =(5o,-i5o,0, , ,0)T, 

V2 go go 

(0,0,ffs,0,0,--)T. 

5s 



^0 - Wnl> 



(6.12) 



Here go 



— (w+ +tj_), 5^ 



= y/muJs ■ The solutions aj(t) of (|4.7I) are found by (|5.8|l and are 



normalized by condition (|4.8|) . They form a matrix A{t) H4.13|l which we rewrite in the block form 



A{t) 



B{t) B*{t) 
C{t) C*{t) 
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(6.13) 



where matrices B{t), C{t) have the foUowing form 

\ 

, C{t) 







V2 










V2 














V 



V^go V^go 




J 



Using 
form 



{ , H5.3()|l , and 14.11|l we find a solution of the last eq uation of system (|4.5|l in the block 
in terms of the blocks B{t),C{t) of the matrix A{t) Ht).13|l as follows 



^ppit) = i^[B{t)D{v)B+{t)+B*{t)D{v)B^t) 
a,,{t) ^^{B{t)D{u)C+{t) + B*{t)D{u)C^t)), 



where the diagonal matrix D{v) = diag(z^i + 1/2, + 1/2, + 1/2), i^i,i^2,i'3 = 0, 00, and the 
symbol (*) denotes the Herniitian conjugate matrix. Matrices Uxx, <^pp are diagonal and their 
explicit form is as follows 



<yxx{t) = — diagf 
m V 



Vl+ V2 + '^ Z/1 + Z/2 + 1 2l^3 + 1 



'^pp(t) — diag(^(cc)+ +uj^){ui + V2 + 1), {^+ + +V2 + 1), 2cJs(2t/3 + l)j. 

The non-zero elements of the matrix <Txp {t) ar e cTpixalO — ~'^P2^i — f^i^i ~ ^2)/2. 

After substitution (|6.12|) into (|5.1U|) . (|5.11|) and taking into account that for V defined by H6.2|l . 
the vector V^i^z^^^zo) — 0, and the only non-zero elements of the 6x6 matrix Vww{zq,zq) are 



= — Vo/7^ for i = 4, 5, 6, we obtain the energy spectrum Ey of the Hamiltonian (ll.l|l . 



UJ+ + LO- 



1 



2ws 



(6.14) 



Note that in the case of zero magnetic field i7 = 0, a similar expression for spectrum was obtained 
in [40]. 

6.2. One-dimensional case 

Consider equation (|l.l|l with linear operator 7i{t) in the form 



" 2to ^ 2 ' 



and the nonlinear operator V(^{t)) as follows 



+00 



V{^{t))^{x, t) = / V{x, y)\^{y, t)\'dy^{x, i), V{x, y) - exp 



{x - yf 
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(6.15) 



(6.16) 



In the absence of a magnetic field {H = 0), the cyclotron frequency (|6.3(l is equal to zero [wh = 0), 
and thus from 1)6. 4|l we find LOa = ^o- Then for the Ritz frequencies we have 



^0 - Wnl- 
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The Ehrenfest system with accuracy 0{T?/^) for operators H6.16|) has the form 



p — ~kx, 

P_ (6.17) 



X = 
m 



'XX — '^xpi 

m 



1 2 (6-18) 

XX 1 



^ xp — ^ pp 

m 

n 2 



As a stationary solution of subsystem H6.17() we take the zero solution 

Zo(n) = {Po{h),X^{h)y = (0,0)T, (6.19) 

In this case matrix A{t) H6.13|l of the pseudo-system-in-variations (|4.7|l is a 2 x 2-niatrix whose 
scalar blocks B{t) and C{t) satisfy equations 

B = -muj^C, C (6.20) 
Floquet solutions ijHT^ a{t) = (^(t), C(t)) ^ of the syste m in variations (|6.12l) normalized by 



condition (a, J'^a*) — 2i, J — i ^ 0^ ) '^^^ written in the form 



ait) 



exp(iLJst) f imujs 
1 



Then for the energy spectrum £"„ of the Hartree-type equation (|1.1() . H6.15|l . (|6.16() we obtain 



Concluding remarks 

An approach to the problem of correspondence between classical and quantum models in the 
nonlinear case has significant differences from the one feasible in linear quantum mechanics. In the 
linear quantum case, a transition from quantum to classical system in the sense of Ehrenfest requires 
a certain property of a quantum-mechanical solution, namely the function has to be trajectory- 
coherent HU.5|I . A state which does not obey this condition is considered to be essentially quantum, 
but one which obeys it is near-classical. For the near-classical states the classical dynamics obtained 
in the limit fi. — > is defined by the classical Hamilton function, and appears to be the same 
regardless whether the quantum solution is localized (at each moment of time) at a point, on 
a curve or on a surface. (Exact meaning of localization on a curve or surface is explained e.g. 
in [2,3].) 

For the Hartree-type equation the situation is different. Classical equations (lU.lOf) (or (|3.9I) ) are 
valid only for states concentrated at a point in each moment of time. These classical equations are 
distinct from those obtained in [55-58] . The latter are integro-differential equations which describe 
dynamics of the n-dimensional manifolds in the 2n dimensional phase space. It was shown that 
for Hartee-type equations an implementation of Horn's approach leads to those integro-differential 
equations for characteristics of a non-local (Vlasov) equation which describes the evolution of the 
classical density matrix. (Recall that in the linear case, the classical density matrix obeys the local 
Liouville equation, whose characteristics are trajectories of classical mechanics.) Thus dynamics 
of point-wise and elongated objects have different equations in the case of Hartree-type models. A 
rigorous derivation of classical equations describing dynamics of fc-dimensional objects (0 < fc < n) 
in 271-dimensional phase space constitutes a separate open problem. 

The Ehrenfest systems (|3.5f) are subjects of mathematical interest, independent of their quan- 
tum origin. Questions similar to those in the linear case, such as about their Poisson structure, 
and the stability of the solutions (including stability with respect to the nonlinearity parameter 
>c) can be addressed in future study. 

In this paper we have shown how the energy spectrum for the Hartree-type equation can be 
retrieved from a rest-point solution of related Ehrenfest system. Similarly, other quantum char- 
acteristics such as quasi-energy spectrum, geometric and adiabatic phases, can be reconstructed 
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from solutions of corresponding Ehrenfest systems. We will attempt to show that in details in our 
future publications. Note that in our approach the quantum characteristics can be found without 
solving the quantum equation. This is particularly valuable and advantageous due to lack of gen- 
eral methods for solving Hartree-type equations. 
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Appendix A 

statement A.l. A solution of equation Hl.l|) with the same Hamiltonian as in (|5.3|l and initial 
condition ^{x,t)\t=f) ~ ijj E V'^{zq) H2.11|I may he written in the form (|5.t)|l . where {</'i/(^)}j^|=o' 

^Pv G T^hi^o) is a complete set of stationary trajectory- concentrated states of equation (|5.3|l . with 
corresponding eigenvalues E^, and S[Q^{t^h)\ is defined in (|5.2|) 

This statement follows from the nonlinear semiclassical superposition principle, which we derive 
first. 

Lemma A.l. Let < "^^(x,t) > be a set of semiclassically-concentrated mod h?''^ solutions of 
equation ()l.l|l with initial conditions 5'i/(a;,0) — il^,j{x) £ 'P^{zo). Then function 

N 

*(x,t) =exp(^5[£l^(t,ri)]) J2 C.exp(-^5[0^Ji,n)])*,(a-,t) + O(?i3/2), (A.l) 

kl=o 

also is a semiclassically-concentrated mod h'^/'^ solution of equation with initial condition 

N 

\E'(a;, 0) = i^{x) — C,jipv[x) G 7'jJ(zo)- Here notation (|5.2f) was used. 

\u\=0 

Proof. Consider associated mod fi.'^/^ to (|1.1|) the linearized in V\ Schrodinger equation [36-39] 

d ^ \ 

-z/i— +.^o(flv(^.^)))* = 0. (A.2) 

Bo(flv(*'^)) = {^{z.w) + ^Sp [K;^(z,u;)A2] + (^),(z),Az) + 

+ i(Az,i5,,(z)Az)| , M = z-Z[t,h). (A.3) 

Z J w—z — Z{t,h) 

Here g^{t, h) = {Z{t, /i), A2) is defined in (|4.2|l . vector i3z(z) and matrix Sjzz{z) are defined in (|3.8|l . 
Substitute for the argument of Sjo{Q^{t, ft)) IIA.3|) the expansion s{t,h)=Q^'^Ht,h)+hQ^^Ht,h) (|OJ, 
where Z^^^^ (t) H4.4|l is the principal term of the phase space trajectory. Then with precision 0{h^^^), 
we obtain 

^o(fl^(i,fiO) +(-^.(^'"H<)),Azo) + i(Azo,i3..(^("Ht))Azo), (A.4) 

where Azq = z — Z^'^'>{t), and io^-* ft)^ is defined by the following equation 

^L'^ (2^{t, ^)) = {^(^, ^) + mU^), Z^'Kt)) + I Sp w)Af\t)]} (A.5) 

A solution mod h'^^'^ of the equation (IA.2|) may be written in the form 

^(^x,t;g^{t,n)) =exp(^5[0v,(t,^)])x(^,i,^'"Hi)), (A.6) 
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and the equation for x^a?, t, , taking into account (|5.2fl . becomes 

- + Axo) + i(Azo,i3..(^(''Ht))Azo))x = 0, (A.7) 

Note that equation l|A.7|) is defined only by trajectory Z'^^\t)^ it is hnear, and it describes evolution 

of any initial state from class 'Pf^{zt^) with = Z'^*'^(0). Thus if the initial state is represented as 
a linear combination 

N 

x{x,Q,zq) = ^ CuXu{x,Q,zn), Xi.(x, 0, zo) € ^"(zo) (A.8) 

\u\=n 

then 

N 



x(x,t,Z(o)(0) = C,X.{x,t,z("\t))+0{h^^'). (A.9) 

|i/|=0 

Now, from (|A.6p we have "^(^x, 0; ^ = x{x, 0, zq), and thus HA.8|I can be written as 



N 



*(f,0)= ^ a*,(f,0), *,(f,0) e7'r°(^o). (A.IO) 

kl=o 



T o com plete the proof o f the lemma it remains to observe that taking into account ljA.6p , equa- 
tion HA.1() is equivalent to (|A.9|I . 

To justify Statement A.l we take N = oo and 

r z 1 
t) = exp - -E^t ip^{x,h). 
I n i 

Then from (|A.1|I we obtain H5.6|l . 

Appendix B 

Here we derive properties of the trajectory-coherent functions, listed in Section 2. 

1. Proof of l|2.2|l . Rewrite the Weyl symbol of the operator {Az}" in the form 

(Az)" = (Ap)"" {Axf^ , (ap, a,) = a. 

Thus, in accordance with H1.5|l . we obtain the following formula for the mean value aa{t, h) of the 
operator {Az}" 

ac.{t,h) - ($|{Ai}"|$) = J dxdydp^*{x,t,h)x 

X expQ(f-y,p))[Apr'(^^y^) ^ mt.h). 

Here 

Ay = y-X{t,h). 

After a change of variables 

Af = Ay = VhC, Ap^Vhuj 
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and using the formula for function t, ft) from the class 'P^{Z{t, h)) H2.1() . we find 

'^"^^'^^ " .^--i^ft3"/2?il"l/22-|-x| J d^dCdiOip*{Ct,h)x 

X exp{z(e- C, + Cr''^{(,t,h) = 

||$||2 = ri"/2 / d^V(f, fiM^ t, h) = n^/^pioit, n). 



Recall that the function if{^, t, K) depends on \fh regularly, and Mo{t, ft) > 0. Therefore 



^ ' Mo(t,ft) te[o,T] Mo(i,ft) ^ 



Q.E.D. 



2. Proof of (|2.3|) follows from the explicit form of a trajectory-coherent function ^{x, t, ft) € 
l|2.1|l and the estimations H2.2|l . 

3. Proof of H2.5|l . Consider a function (/)(x) G §. Then for any function ^{x,t,h) E V\ the 
integral 



dx 



after the change of variables ^ = Ax / Vft becomes 



(|a>(t,ft)|2|^) = — / 0(x(i,ft) + yftO|^(c,i,ft)|^de 



|^(t,ft)| 

Taking the limit ft — > and using that 

||^(t,ft)|P=ft"/2 / |(^(C,<,ft)|2de 



where the function (p[$^,t,h) depends on \/H regularly, we obtain the statement. 

Proof of H2.6|l is similar to the previous one, if we note that the Fourier image of the function 
^{x, t, ft) G may be represented in the form 

t, ft) = exp{^ [^(t, ft) - (p, X{t, ft))] }^{ ^^^'^\ t, hj , 
where 

mt,h)^^^Je-^(^'^^^{U,h)dC 
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